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Abstract 

We have analyzed the B n and Bg ^ K semileptonic form factors and B — > V'j{V = 
K*,p,uj) processes in the hght-cone QCD sum rules. In order to enhance the predictivity 
and reliability of numerical results the chiral-current correlator is employed and the twist-3 
light-cone wavefunction can be effectively eliminated from the sum rules. 
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1. Introduction 



B decays into a light meson are an important ground to understand and test the standard 
model(SM), since they can provide the signal of CP- violation phenomena and perhaps, a window 
into new physic beyond the SM. Nevertheless, the definite interpretation for the relevant ex- 
perimental data demand that we have the ability to precisely compute the physical amplitudes. 
One is forced to use some approximate methods, such as QCD factorization formula fll], heavy 
quark effective theory(HQET), chiral perturbative theory(CPT), Modified PQCD approach[0, 
QCD sum rules@] and light cone QCD sum rules0, ||. Each of them has advantages and disad- 
vantages. QCD factorization formula is viewed as a great progress in phenomenology of heavy 
fiavors, however the underlying long distance effects included in a series of hadronic matrix 
elements confront still us, which are not being dealt rigorously with and thus would hinder 
us from doing such a desired calculation. CHPT and HQET, as two effective theories at low 
energy, can describe very well light-to-light and heavy-to-heavy transitions, respectively. They 
are not suitable for a study on heavy-to-light processes. It seems that the modified PQCD 
approach with a resummation of Sudakov logarithms is consistent with the physical picture of 
heavy-to-light transitions due to the hard gluon exchange. However, a detailed analysis shows 
that the reliable PQCD calculation depends on whether the singularities can be eliminated or 
suppressed. The singularities come from on shell gluon, on-shell light quark and on-shell heavy 
quark. QCD sum rule method is built on the basis of QCD theory, but the obtained form 
factors behave very badly in the heavy quark limit Mq — > oo due to the fact that one omits 
the effects of the finite correlation length between the quarks and the gluons in the physical 
vacuum. In order to overcome the defect, QCD light cone sum rule(LCSR) approach has been 
developed in Ref. 0, |] and has been regarded as an advanced tool to deal with heavy-to-light 
exclusive processes. In comparison with the case of the traditional QCD sum rules, in LCSR 
approach the operator product expansion(OPE) is carried out near the light cone x"^ ~ 0, in- 
stead of in small distance a; ~ 0, and the non-perturbative dynamics is parametrized in term of 
so-called the light cone wavefunctions of light mesons, instead of the vacuum condensates. It 
is the striking advantage of the approach to describe heavy-to-light transitions in a way consis- 
tent with the universally accepted viewpoint that non-perturbative QCD dynamics occupies an 
dominant place and perturbative hard gluon exchanges contribute only a subleading effect in 
that case. Very recently, it was generalized to study the non-factorizable effects in i? — >■ vrvrP 
and to probe heavy-to-light form factors in the whole kinematically accessible ranges]^. The 
technical details of LCSR can be found, for instance , in ||^, while for a detailed comparison 
with the traditional sum rules, see [^]. However, a problem with the approach is that the rel- 
evant light cone wavefunctions are intrinsically non-perturbative and can only be determined 
via some phenomenological methods. This could considerably affect the reliability of sum rule 
predictions. Accordingly, it is essentially important to find out some certain way to control the 
pollution by the light cone wavefunctions. With this end in view, a pragmatic strategy has been 
suggested in Ref. [5,10-13], in which some specific chiral current operators act as the interpolat- 
ing fields in the correlators used as a LCSR calculation on the form factors for heavy-to-light 
transitions, making some twist-3 components cancel out in the OPE. Taking it into account 
that these twist-3 distributions are equally important in most cases but poorly known to us 
compared with the corresponding twist-2 ones, we may think of this procedure as an effective 
way to narrow down the uncertainties in LCSR calculations. Here we intend to present a simple 
summary of the our previous works [10-13] on applications of the improved LCSR approach to 
some important heavy-to-light transition processes. 

2. Heavy-to-light Semileptonic Form Factors [|TD|, |TTj] 

Semileptonic B decays included by the h ^ u transition are regarded as the most promising 
processes adequate to extract \Vub\ from the relevant data. We will focus on a discussion on 
B ^ TT and B ^ K semileptonic transition form factors within the framework of an improved 
LCSR.. 



The form factors for S — > tt transitions can be defined as follows, 

{n{p)\u^,b\B{p + q)) = 2f--{q^)p, + f^--{q')q„ (1) 
with q being the momentum transfer. We choose to use a chiral current correlator 

n^(p,?) =^ / ciW''"(7r(p)|T{M(a;)7^(l + 75)6(x),6(0)z(l + 75)d(0)}|0) 

= ^(?^ (p + qf)p, + U{q', {p + qf)q,, (2) 

to calculate f^^^{q^) and f^^^{q^)- Its hadronic expansion reads 
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We argue that such a procedure is conceptually reasonable and technically feasible by ob- 
serving the hadronic form of (2). Because of the special chiral constructs , the correlator would 
receive the contribution of the 0"^S mesons, in addition to the Q~B meson, in its hadronic 
expression. However, this causes no difficulty in extracting the desired sum rule for the form 
factors. Since the lowest Q'^B meson is far from the ground state of Q~B mesons and slightly 
below the first excited Q~ B mesons in mass so that we can safely isolate the pole contribution 
of the lowest Q~B meson and parametrize them from the higher 0~ states as well as all the Q^B 
mesons in a dispersion integral: 

(p + qf) = ^ J pH^) ^ subtractions, (4) 



m6(m| - (p + J s-ip + qY 



so 



n^(g', (P + qf) = /^f^'^Bh _^ 7 P"{s) , ^ subtractions. (5) 
mb{m% - {P + q)) J s-{p + qy 

Here the threshold parameter Sq should be set near the squared mass of the lowest O'^B meson. 

QCD calculation of the underlying correlator may be allowed, on the other side, for the 
negative and large p^ and (p + q)^, in which case the OPE goes effectively in powers of the 
deviation from the light cone x'^ Ri 0. The chiral hmit p^ = = will be taken throughout 
this discussion, for simplicity. Carrying out the OPE for the correlator, we have 



YlQCD^^2^ (p + q)^) = H('^«)[g^ (p + qf] + U^'i'^^^q^, (p + q) 
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with n^^') being the two-particle contribution and n'-''''^) the three-particle one. (fniu) is the 
leading twist-2 wavefunction and the others have twist-4. They are defined as follows 
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At this point, the important observation is that difference from the standard LCSR calcu- 
lations, in the present case the twist-3 components cancel precisely out in the OPE. It can 
improve greatly the precision of LCSR predictions. 

Now the LCSR for /(g^) can be obtained using the standard procedure. The result is 
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To further proceed, we need to make a choice of input parameters. The parameters entering 
the sum rule are the h quark mass m^, B meson mass m^, decay constant and threshold 
parameter sq. We take = 4.7 — 4.9 GeV, = 5.279 GeV. As for the decay constant fs 
and the threshold parameter sq, the two-point correlator 

K{q')=t [ ciW«-(0|g(x)(l + 75)Ka:),&(0)(l-75)g(0)|0) 



will be used to estimate it, for consistency. A standard manipulation yields three self-consistent 
sets of SR results: (1) fe = 165 MeV for = 4.7 GeV and Sq = 33 GeV^ (2) fs = 120 MeV 
for mb = 4.8 GeV and Sq = 32 GeV^ and (3) fe = 85 MeV for m^ = 4.9 GeV and 
So = 30 GeV^. The parameters relevant to the vr meson contain the decay constant and the 
set of light cone wavefunctions. We use = 0.132 GeV and the wavefunction models |T^, |TB 
based on the conformal invariance of QCD: 



3 3 



with the coefficients a2{u.Ub) = 0.35 and 04 = 0.18, at the scale Ub = (m^ — mlY^'^ 
2.5 GeV, and the Gegenbar polynomials C2^^{2u - 1) = f [5(2m - 1)^] and Cf ^(2m - 1^ 
[21(2u - 1)^ - 14(2m - 1)2 + 1], and 
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+10?f lnM(2 - 3m+ -m^)], 
o 

g2{u, Ub) = ^6'^uu{u - u). 

^^{ai) = 306\ai - a2)al^ + 2e(l - 2a3)], 

ip^iai) = m^alil - a3)[^ + 2e(l - 2^3)], 
i^wiai) = 1205'^e{ai - Q;2)aia2a3, 

= -1206^aia2a3[^ + e(l - 3^3)], 

with 5^{iXb) = 0.17 GeV^ and = 0.36. 

Having fixed the input parameters, we can carry out the numerical analysis. The reasonable 
range of is found to be 8 GeV"^ < < 17 GeV"^, in which the variation of /(g^) with 
turns out to be negligibly small. For a specific = 12 GeV^, the sum rules for f^^'^{q^ = 0) 
are predicted to /(O) = 0.27, 0.29 and 0.33 (corresponding to the set (3), set (2) and set 



(1), respectively), which are in basic agreement with an estimate from the standard LCSR|H|. 
As a matter of fact, numerical agreement between the two different approaches exists up to 
= 10 GeV^ the differences being within 20%; while for 10 GeV^ < < 18 GeV'^ there is 
a numerical disagreement of about 20 — 25%. For the region > 18 GeV^, application of the 
LCSR is questionable such that a comparison is meaningless between the two sum rule results. 

The total uncertainty in f^^^iq"^), which arises from the uncertainties in the b quark mass 
nib, the threshold parameter sq, the decay constant fs and the pionic light cone wavefunctions, 
is estimated at the level of about 26%. 

Following the same procedure, we may analyze the Bg ^ K form factor /^"^^(g^). As 
compared with the case of B —>■ tt, however, the f^^^^iq^) is more difficult to calculate, for 
SU(3) breaking corrections to the twist-3 wavefunctions of K meson have not been investigated 
completely in the literature. Explicitly, this problem can be avoided in our approach. 

The correlator used for the LCSR calculation on f^^^^{q^) may be obtained by an obvious 
replacement d ^ s m Eq. (2). Carrying out the OPE and considering the non- negligible 
i^-meson mass we have 

n"tf . (p + ,f) = VEirWsjB. l^piL-is (12) 

^ ^ {nib + ms){m\ - {p + qf) J s - {p + q^ ^ ^ 

" so 

where the definitions P = ai + aa^, A = {mf — q^)/{so — q^ — ni\) and Dai = daida2da36{l — 
a\ — a2 — cx-i) have been used, and the light cone wavefunctions of the K meson obey the same 
definitions as the corresponding cases of tt meson. 



In the following numerical analysis we adopt a model presented in 



<^k{u) = 6u{l - m) {1 + 1.8 [{2u - 1)2 - 1/5] - 0.5(2n - 1) [l + 1.2[(2n - 1)^ - 3/7]] } (13) 



for the leading twist-2 wavefunctions of the K meson, and neglect the SU(3) breaking effects 
for all the twist-4 distributions. Numerically, we take fx = 0.16 GeV and = 0.15 GeV. 



Furthermore we use mb = 4.8 GeV, ttib^ = 5.369 GeV, rriB* = 5.325 GeV, sq = 34 GeV , 
/b, = 0.142 GeV, and /b* = 0.132 GeV in the B channel. 

With these inputs, the fiducial interval of is determined to be 8 < < 17 GeV^, 
depending slightly on g^, for = — 17 GeV^. The sum rule results for /^"^^{q'^) are 
illustrated by the solid line in Fig. 1, the resulting total uncertainty being about 20%. 

Having in mind that all the above discussions are limited to the region of low and intermediate 
momentum transfer, for which the OPE is effective, one must find another way to estimate the 
form factors at the large momentum transfer. The heavy-to-light transition form factors in 
the whole kinematical accessible range, for example, 
precisely be represented as 



Bs K form factor f^''~'^{q^) can be 
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2/77,^* (1 — q^ /m\ 
= FG{q') + FHiq'), 

with gB*BsK being the strong coupling defined by {B*{q, e)K{p)\Bs{p+q)) = —gB*BsK{p-e), and 
p(cr) a spectral function with the threshold (Tq. Obviously, Fci^q^) stands for the contribution 
from the B* pole, which describes the principal behavior of the form factor around q^ = q'^^x^ 
and Fui^q^) parametrizes the higher state effects in the B* channel. The non-perturbative 
parameter fB*gB*BsK is calculable with the correlator used for the LCSR estimate of fic^^i'f)- 
Accordingly, modelling the higher state contributions by a certain assumption and then fitting 
Eq.(14) to its LCSR result f^^^^{q^) in the region accessible to the light cone OPE, we might 
derive the form factor f^'^^iq"^) in the total kinematical range with a better accuracy. 

We work in the large space-like momentum regions (g^ <^ and {p + g)^ <^ 0) for the 
correlator in question. Making the Borel improvements on the yielded theoretical expression 
F'^^^(q'^, (p + qY) — > F'^^^(Mf, M2), and then matching it onto the corresponding Borel 
improved hadronic form via the use of quark-hadronic duality ansatz, the final sum rule for 
fB*9B*BsK reads. 
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(15) 



The resulting sum rule for fB*gB*BsK is numerically fB*gB*BsK = 3.88 GeV, with the uncer- 
tainty of about 18%. Using its central value, we obtain a B* pole approximation for /'^^^^(g^), 
which is shown by the dotted line in Fig.l. It is explicitly demonstrate that a perfect match 
between the direct LCSR and B* pole predictions appears at g^ 15 — 20 GeV^. 

It is important and interesting to make a comparison of our sum rule results and those from 
the standard LCSR based on the correlator vector and pseudoscalor currents, which are easy 
to obtain using the twist-3 wavefunction suggested in Ref. |T^, leaving the twist-4 distribution 
amplitudes unchanged and making a corresponding replacement of the other relevant input 
parameters in Eq.(79) and (44) of Ref. [Q.We observe that the standard approach gives the 
same matching range as in our case and the resulting deviations from our predictions turn out 
to be between —10% 15%, depending on g^, in the total kinematically accessible region. 

Assuming the higher state contribution in Eq.(14) to obey Fui^q^) = a/ (1 — hq^ /m\^ — cq^/m^*), 
we could give a model for the form factor 



fB*gB*BsK 
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,4 I. 



(16) 



The parameter a can easily be fixed at —0.07, using the central values of filT^^i^) 
fB*gB*BsK- In the region = — 18 GeV"^, the best fit of Eq.(16) to fLc{(f) yields h = 1.11 
and c = —8.33. It turns out that the fitting results(the dashed line in Fig.l) reproduce precisely 
the LCSR prediction up to = 18 GeV^ and support considerably the B* pole description of 
the Bs K form factor in large region. 

Finally, we look roughly into SU(3) breaking effects in heayy-to-light decays by considering 
the ratio of the derived f^"~*^{q'^) over f^'^^iq"^), which can be modelled using all the same 
method as in the Bg K case. For the common kinematical region the resulting ratios, 
a comparable result 1.05 — 1.15 with that from the standard approach, favor a small SU(3) 
breaking effect. 

3. The rare decays B {K*, p, cu) + 7[]g, 



The rare decays B {K*, p, u) +7 induced by the flavor- changing neutral currents (FCNC's) 
provide us with a good opportunity to search for new physics. These decay modes are dominated 
by the electromagnetic penguin operators responsible for b (s, d) + 7. The relevant form 
factors F^^(-^*'^'''^)(0) have been studied widely by using the various tricks. With the standard 
LCSR, ones found F^^^' (0) = 0.32 + 0.05 ^ and ^^^''(O) = 0.285 ± 15%[|l|. 

Here we reanalyze them using the improved LCSR approach, to eliminate the pollution by 
some of nonleading distributions. First, we focus on the case of i? — > K*^. The relevant decay 
amplitude reads 

A{B ^ K*^) = Cm^e^ < K*{p,r])\sa^,{l + ^,)q%\B{p + q) >, (17) 

where e and q are the emitted photon polarization vector and momentum, respectively. The 
constant C depends on the CKM matrix elements Vf*Vtb and its apparent form can be found 
in Ref.g. 

The hadronic matrix element in may be parametrized in terms of the form factor F^~*^* (0), 
< K*{p, 7])\sa,,{l + l5)q''b\B{p + q) >= (-2^e^,„^r/*'^gV + 2p ■ gr/^ - 2g ■ r/*p^)F^^^*(0).(18) 
A correlator, which is suitable for our purpose, is 
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The hadronic expression for the invariant function F [{p + g)^] is of the following form 
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so 



Here the dispersion integral stands for the contribution of both higher pseudoscalar states B^ 
and scalar resonance states B^. 

On the other hand, for the region of large spacelike momenta (p + g)^ ^ the resulting QCD 
form of F [{p + g)^], to the twist-2 accuracy, reads. 



U 



s — {p + qY u'^mj^, + m\ 



:ds. 



(21) 



The variable m, the fraction of the K* meson momentum carried by the s quark, is related to 
s by s = mllu — (1 — u)m\*. The leading twist-2 wavefunction ip^iu^uD parametrizes the 
nonlocal matrix element < K*{j),vi)\s{x)a^yq"{l + 75)^(0) |0) > as the following, 



<ir*(p,r/)|s(x)a^,g'^(l + 75)^(0)10 > 



K* 



X / rfMe™f^(/?x(M,M^), 



(22) 



with the decay constant f^* = 210 MeV. 
The final LCSR for the form factor /^*(0) is 
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(23) 



which, in comparison with that in Ref.[|l^], does not receive the contribution of distribution 
amphtudes ipp and g^, and thus is more rehab le. 
Using as inputs the Chernyak-Zhitnitsky(CZ) model for ip^ {u,fi'^), 



(f^ (u, /i^) = 6m(1 — u) 



1 + 0.57cu - 1.35(c<j2 OAGilu^ - u] 



(24) 



with u = 2u — l, and taking into account all the possible uncertainties in the numerical analysis, 
we have F^^^*(0) = 0.34 ± 0.05, slightly greater than that in Ref.g and B{B K*-f) = 
(5.1 ± 1.7) X 10~^, a result comparable with the experimental observation pl| Br{B K*'~f) = 
(4.5 ± 1.5 ±0.9) X 10-5. 

The corresponding relations in i? ^ (p, cu) +7 cases can be obtained by making replacements 
s ^ d, ip^* {u, jj,'^) — i> Lp^l''^\u, jj?), rriK* ^{p,lo) and f^* f'f''^^- The isospin symmetry 
allows us to adopt the same wavefunction model for p+, p", p° and ui. We choose the two 
different models 



V?±(m, p^) = 6m(1 — u) 



1 + 0.138 X 75(cj2-0.2) 



and 

V9j_(u, p^) = 6m(1 — u) 



1 + 0.077 X 7.5(0;^ - 0.2) - 0.077 x (39.375a;^ - 26.25a;^ + 1.875 



(25) 



,(26) 



which are respectively introduced by Ball and Braun in Ref. p2[ and Bakulev and Mikhailo 
in Ref. |23|, and called the model I and the model II. In the two cases, the resulting sum rule 
predictions are summarized as follows: 



and 



F^-^^*(0) = 0.335 ± 0.050, Br{B ^ p±7) = (2.71 ± 1.00) x 10"^ 

(p°, cu) + 7) = (1.36 ± 0.50) X 10"^ 

p±7) = (1.79 ±0.61) X 10"^ 



F^^^°/'^(0) = 0.237 ±0.035, Br{B 
for the model I, and 

^B-P±(0) = 0.272 ± 0.029, Br{B 

and 

F^^^°/'^(0) = 0.192 ± 0.021, Br{B (p°, cj) ± 7) = (0.90 ± 0.31) x 10"^ 

for the model II. The results with the model I are slightly great than those in Ref. |jl9|, while 
the ones with the model II accord with those in Ref. ||19||. Compared with the experimental 



(27) 
(28) 

(29) 
(30) 



observation Ref.H Br{B^ PO7) < 3.9 x 10-^ Br{B^ CU7) < 1.3 x 10"^ and Br{B- 
P^l) < 1.1 X 10~5, our predictions are below the experimental upper limits. 

4. Conclusion 

Light-cone QCD sum rules is a good framework for calculating heavy-to-light form factors. We 
have analyzed the. B —>■ tt and Bg K semileptonic form factors, also B —>■ K*^ and 
B (p, C(j)7 processes. In order to enhance the predictivity and reliability of numerical results 
it is worthwhile to study how to reduce the various uncertainties in the light-cone QCD sum 
rules. Here we employ the chiral- current correlator. It is explicitly shown that the twist-3 
light-cone wavefunction, which have not been understood well, can be effectively eliminated 
from the sum rules. 
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Figure 1: The Bg ^ K form factor /{q^) in the total kinematical range. The sohd hne denotes 
the LCSR result fLc{(f), which is rehable for < < 17 GeV"^. The dotted hne expresses 
the B* pole prediction suitable for large q^. The best fit of Eq. (11) to fhciq^) is illustrated 
by the dashed line. It should be understood that the plotted curves correspond to the central 
values of all the relevant parameters. 




